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We calculate pZ^ potentials from the equal-time Bethe-Salpeter amplitude measured in the 
quenched QCD simulation with the spatial lattice volume, (4.4 fm)^. The standard Wilson gauge 
action with the gauge coupling j3 = 5.7 on 32"* lattice together with the standard Wilson quark 
action are used. The hopping parameter K^d — 0.1678 is chosen for u and d quarks, which cor- 
responds to m-,1 ~ 0.37 GeV. The physical strange quark mass is used by taking the parameter 
Ks = 0.1643 which is deduced from the physical K meson mass. The lattice spacing a = 0.1420 
fm is determined by the physical p meson mass. We find that the pZ^ potential has strong spin 
dependence. Strong repulsive core is found in '^o channel while the effective central potential in 
the ^Si channel has relatively weak repulsive core. The potentials also have weak attractive parts 
in the medium to long distance region (0.6 fm < r < 1 .2 fm) in both of the '^o and ^Si channels. 
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1. Introduction 

Study of hypernuclei is one of the frontiers in nuclear physics. The strange degree of freedom 
gives a new dimension to the description of nuclear structure. A hyperon (or a strange quark) em- 
bedded in nuclei plays a characteristic role as an "impurity" or a probe in the many body system[|l|]. 

Modern nucleon-nucleon (A^A^) potentials give a successful description of the A^A'^ scattering 
data and have been used to make precise calculations of light nuclei |§, ^. In contrast, hyperon- 
nucleon {YN) and hyperon-hyperon {YY) interactions have large uncertainties, because the scat- 
tering experiments are either difficult or impossible due to the short life-time of hyperons. The 
full phase shift analysis at the same level of the A^A'^ case is not available yet. Although a lot of 
theoretical models which describe FA'^ interaction together with the A^A'^ interaction have been pub- 
lished, the different model predicts different phase shifts and scattering parameters, e.g., of the AA'^ 
potential in spite of the nice description of the A^A'^ sector [^, ^, ^, ||, ?]. 

Recently, a lattice QCD study of the A^A'^ potential has been performed[p^]. This approach may 
lead to a new paradigm to study the FA'^ and YY interactions too, since the lattice QCD is an ab 
initio method of treating the fundamental theory of strong interaction. (See also the conventional 
approach to the FA^ phase shifts using the Luscher's finite volume formula[jn|].) 

The purpose of the present report is to explore the YN and YY potentials from the lattice QCD 
simulation on the basis of the methodology developed in Refs. 1 1^, 12]. The extension from the A^A'^ 



potential to FA'^ and YY potentials is relatively straightforward. In the case of the A^A'^ potential, there 
are only two representations in the isospin channel, i.e., 2(E)2 = 3©lin SU (2), which correspond 
to isovector and isoscalar channels. Including the strange degrees of freedom extends the arithmetic 
into flavor (3), 808 = 27 010010801008. Here the isovector (isoscalar) channel of the A^A^ 
sector is assigned to be a subset in the 27-plet (10-plet) representation. The potentials for newly 
arising channels are hardly determined from the real experiment so far. The lattice QCD simulation 
with physical strange quark mass provides new numerical "data" for these strange channels. 

In this report, we focus on the A'^S potential in the isovector (7=1) channel as a first step. 
We keep away, at present, from the isoscalar (7 = 0) channel of the A'^E potential, since it is not 
the lowest state of the isoscalar '^o channel and AA strong decay mode may open below the A^S. 
There are almost no experimental information on the A^S interaction, although a few experimental 



data[]13|, [14|, |lj] suggest that the S-nucleus potential would be weakly attractive. Moreover, the 
S-nucleus interaction will be studied as a day one experiment in the near future at J-R\RC[p^ 
through {K^,K^) reaction on the nuclear target such as ^^C. 

2. Formulation 



As we mentioned, the methodology to obtain the potential is along the lines of Refs. [12] 



and [10|. The latter describes, successfully for the first time, the A^A'^ potentials from lattice QCD 



simulation. We start from an effective Schrodinger equation for A^S system at low energies: 

-2^V20(r)+ J d'r'U{f,l')^{f')=E^{r), (2.1) 

where jJ. = 171^111^/ (m/v + m^) and E are the reduced mass of the A'^S system and the nonrelativistic 
energy in the center-of-mass frame, respectively. The nonlocal potential can be represented by the 
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derivative expansion as 

U{f,f)=VN^{r,W)5{f-r). (2.2) 
The general expression of the potential V^z would be 

Vnz = Vo(r) + Va(r)(a^ • d.)+V,{r){% ■ T^) +y^,(r)(a^ • d.){TN ■ ^s) 

+VT{r)Sn + VT,{r)Sn{'CN ■ r.) + VLsir)iL ■ S+) +VLSrir)iL ■ S+){Tn ■ Tz) 
+VALs{r){L ■ S-)+VALST{r){L ■ S-){Tn ■ fs) + ^(V^). (2.3) 

Here Sn = 3(ai • n){d2 - n) — Oi ■ $2 is the tensor operator with n = r/\r\, S± = {Oi ± 02) /2 sym- 
metric (+) and antisymmetric (— ) spin operators, L = —if x V the relative angular momentum 
operator, and % (fs) is isospin operator for N (S). We note that the antisymmetric spin-orbit 
forces (Vals and Valsx) newly come up because the constituents (A'^ and S) are not identical. 

According to the above expansion of the potential, the wave function should be classified by 
the total isospin /, the total angular momentum and parity with J = L + S+. Particular spin 
(isospin) projection can be made in terms of dj^ ■ dz (% • ^s). e-g-» for the isospin projection we 
have = (1 — T/v fs)/4 and p(^=^) = (3 + % • fs;)/4. In this work, we focus only the isospin 

7 = 1, 5- wave component of the wave function, so as to obtain the (effective) central potential 
through 

The S'-wave wave function is measured from the equal-time Bethe-Salpeter (BS) amplitude 
<p{r;k) as 



P^^;k), (2.5) 



Pa{x) = eahc{ua{x)Cysdh{x))uca{x), (2.6) 
sJCj) = £ahc{ua{y)Cy5Sb{y))s,p{y), (2.7) 

where a and jS denote the Dirac indices, a, b and c the color indices, and C = 7472 the charge 
conjugation matrix. The summation over ^ € is taken for cubic transformation group to project 
onto the S-wave, and the summation over x for zero total spatial momentum, pa {x) and {y) are 
the local field operators for the proton and "EP. We take the upper components of the Dirac indices 
a and j8 to construct the spin singlet (triplet) channel by P^^ = {02)ap iP^p = {(^i)aii)- 

The (/>(r) with f = x — y is understood as the probability amplitude to find "nucleon-like" 
three quarks located at point x and "S-like" three quarks located at point y. The <p (r) includes not 
only the elastic amplitude NE — s- NE but also the inelastic amplitudes such as NE — > nNE and 
NE AT, and so on. Note that, at low energies below the thresholds, however, the asymptotic 
behavior of ^(r) is not affected by the inelastic contributions, since they decrease exponentially in 
the asymptotic region. (In the present calculation with the isospin 7 = 1, the AA channel is closed 
due to the isospin conservation.) On the other hand, ^{r) and hence the potential may depend on 



the interpolating fields in the interaction region. Further study on this issue is found in Ref.[17] for 
the NN potential. 
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In the actual simulations, the BS amplitude is obtained through the four-point correlator, 







n ) e 



-E„(t-to) 



(2.8) 
(2.9) 



Here ^^^o (?o) is a source term located utt = to. We utilize the wall source for ^ pi^o in this work 
in order to enhance the lowest scattering state of the pS" system. En is the energy of the pS" state, 
\n), andA„(fo) = («l7^pso(^o)|0). 

We also calculate the two-point correlator, C{t;to) = llx{'^\^ai^,t)^a{x' ,to)\0), for the octet 
baryons = A'^,S,A,r), in order to check whether various two baryon (AA,A'^E,AZ and ££) 
thresholds are reproduced in the correct order. The interpolating fields for A and £+, employed in 
this work, are given by 

Aaix) = eabc{{daix)CY5Sb{x)) Uca{x) + {Sa{x)CysUb{x)) dca{x) - 2 {Ua{x)Cy5db{x)) Sca{x)} , 

(2.10) 

= -£abc{ua{y)CY5Sb{y))ucfi{y). (2.11) 



3. Numerical calculation 

We use the standard Wilson gauge action at the gauge coupling j8 = 5.7 on the 32^ x 32 



lattice together with the standard Wilson quark action. See Ref.[]17[] for details. The hopping 
parameter of Kud = 0.1678 is chosen for the u and d quarks, which corresponds to m^^ ~ 0.37 GeV, 
tUp ~0.81 GeV,andm/v~ 1.16GeV.In order to determine the parameter for the strange quark mass 
(JCs), we first measure the correlators of pseudo scalar and vector mesons by using the interpolating 
fields given by 

-^psW = ^i.co;W75,ai3'?2,cj3W for pscudo Scalar meson, (3.1) 
-^v./cW = ^i,caW7*:,a/3'?2,c/3W for vcctor meson, (3.2) 

with applying six set of hopping parameters; first three sets for Ki = K2 (7t and p mesons) with tak- 
ing a number from {0.1678,0.1665,0.1640}, and another three sets for Ki > K2 {K and K* mesons) 
with taking different two numbers from {0.1678,0.1665,0.1640}. Assuming the following func- 
tional forms for the pseudo scalar meson mass squared and for the vector meson mass, 

2 B f I 1 \ B / 1 1 



D f I 1 \ D f I 1 \ 

we obtained critical hopping parameter Kc = 0.1693, and the physical parameter, Kphys = 0.1691, 
from {m^a/mpa) = (135/770). The lattice scale is determined so as to be a = 0.1420 fm from the 
physical p meson mass. The parameter for the strange quark mass is determined as = 0.1643 
from the physical K meson mass (494 MeV). 
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Figure 1: The radial wave function of pS", in ^So (circle) and (triangle) channels, obtained from the 
lattice QCD at f - fo = 6. 



4. Results and Discussion 

Table |l| lists the hadron masses measured from the present lattice QCD simulation. 1283 gauge 
configurations are used to calculate the hadron masses and wave functions. (17 exceptional config- 
urations are not used out of totally 1300 gauge configurations.) We note that the present results for 
the baryon masses provide the correct order of particular threshold energy of two baryon states in 
the strangeness S = —2 sector; £',/,(AA) = 2525(1 l)MeV (this channel is not allowed in the present 
case because of isospin conservation), Eth{NE) = 2544(12)MeV, EthiAL) = 2575(1 l)MeV, and 
£'f/j(EE) = 2624(1 l)MeV. This warrants the desirable asymptotic behavior of the wave function. 

Figure |l] shows the wave function obtained at the time slice t — to = 6. The '^o (^^i) channel 
is plotted by circles (triangles), which are normalized at the spatial boundary r = (32/2,0,0). All 
the data are taken into account for r < 0.7 fm, while only the data on the x-, y-, and z-axis and 
their nearest neighbors are used to plot for the outer region. As seen in the Figure, the wave 
functions are suppressed in the short distance region, and a slight enhancement is found in the 
medium range region for both ^Sq and ^Si channels. There is sizable difference between the '^o and 
^^i, particularly of suppression in the short distance, suggesting that the channel has stronger 
repulsive core. 



ruji nip nifc m/c* mp m^o mi^+ 

367(1) 811(4) 552.6(5) 882(2) 1164(7) 1379(6) 1263(6) 1312(6) 

Table 1: Hadron masses, given in units of MeV, measured from the lattice QCD simulation. The number in 
the parenthesis shows the errorbar in the last digit. 
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Figure 2: The effective central potential for pZ'^, in the '5o (circle) and ^Si (triangle), obtained from the 
wave function at time slice t — to = 6. The inset shows its enlargement. 



Figure § shows the (effective) central potentials for /^S*' in the ^Sq and ^Si channels. These 
results are still preliminary, since the potentials are obtained by assuming E = 0: The energy should 
be determined by fitting the asymptotic behavior of the wave function with the use of the Green's 



function G{f,E) which is a solution of the Helmholtz equation on the lattice|ll2|, [18|]. Preliminary 
calculation suggests that E would be small negative values for both and channels, similar 
to the case of the A^A^ potential[|l^, |19|]. 

In order to see the ground state saturation of the present results, we plot, in Fig. ^, the time-slice 
dependence of the potential in both of the ^Sq (left-hand-side) and (right-hand-side) channels 
at several radial distances; r = 0.14, 0.20, 0.71, 1.42, and 2.27 fm. We can see that the saturation 
is achieved for ? — fo > 6 within errors. 

The present work is a first step toward the YN and YY potentials from the lattice QCD simula- 
tion. Systematic studies of the various channels such as AN, LN, AA, and so on are all interesting 
and important because they are intimately related not only to the structure of hypernuclei but also 
to the internal structure of neutron stars. We will present such studies in the near future. 
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Figure 3: The time-slice dependence of the potential in the (left-hand-side) and (right-hand-side) 
channel for several radial distance r. 
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